Magnetoelastic Coupling in the Spin-Dimer System TlCuCls 
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We present high-resolution measurements of the thermal expansion and the magnetostriction of 
TlCuCl3 which shows field- induced antiferromagnetic order. We find pronounced anomalies in the 
field and temperature dependence of different directions of the lattice signaling a large magnetoelas- 
tic coupling. The phase boundary is extremely sensitive to pressure, e.g. the transition field would 
change by about ±185 %/GPa under uniaxial pressure applied along certain directions. This drastic 
effect can unambiguously be traced back to changes of the intradimer coupling under uniaxial pres- 
sure. The interdimer couplings remain essentially unchanged under pressure, but strongly change 
when Tl is replaced by K. 



One of the most simple quantum spin system is a 
spin- 1/2 dimer. If such dimers are weakly coupled to 
each other, very rich and fascinating physical properties 
are predicted for various theoretical models and can be 
observed experimentally in suitable materials. For ex- 
ample, the two-dimensional Shastry-Sutherland model Q 
is realized experimentally by SrCu2(B03)2 and its low- 
temperature magnetization as a function of magnetic 
field shows distinct plateaus at certain fractional values of 
the saturation magnetization 0] ■ Magnetization plateaus 
are also observed in the three-dimensional spin-dimer sys- 
tem NH4CUCI3 |3- Such plateaus are, however, absent in 
its iso-structural (at 300 K) counterparts -RCUCI3 with 
R = Tl and K, which both have a non-magnetic ground 
state up to a certain magnetic field . Above this field 
a Neel order with staggered magnetization perpendicular 
to the applied field occurs 0, @ and it has been proposed 
that this transition should be viewed as a Bose-Einstein 
condensation (BEC) of magnons fh. lil ioL IToj . Accord- 
ing to a recent neutron scattering study, the different 
behavior of NH4CUCI3 is connected with two structural 
phase transitions in that compound . Despite their 
qualitative similarity, the magnetic systems of TICUCI3 
and KCUCI3 show pronounced quantitative differences. 
The triplet excitations of TICUCI3 are strongly disper- 
sive, whereas those of KCUCI3 have a weak disper- 
sion 0, 0, 0] . Consequently, the minimum gap A ~ 
8 K is significantly smaller than the intradimer coupling 
J ~ 64 K for TICUCI3, whereas this difference is much 
weaker for KCuCl 3 (A ~ 30 K and J ~ 50 K) USEE!. 
The very different behavior of the -RCUCI3 series shows 
that small structural differences strongly influence the 
magnetic subsystem. Evidence for a strong magnetoelas- 
tic coupling in TICUCI3 is also found in ultrasound and 
NMR data, which indicate that the phase transition of 
TlCuCla has a significant contribution of first-order char- 
acter [l4j, |l5j. A BEC is expected to be of second order, 
but spin-phonon coupling can drive a continuous tran- 
sition into a first-order phase transition [16j . Moreover, 



hydrostatic pressure of less than 0.5 GPa is already suffi- 
cient to close A and to induce antiferromagnetic order in 
TICUCI3 without magnetic field 0, 0] . A microscopic 
understanding of the relevant changes under pressure is 
still missing. In particular, it is not clear why external 
pressure decreases A, whereas the substitution of the Tl + 
ions by the smaller K + ions increases A. 

We present high-resolution measurements of the ther- 
mal expansion and the magnetostriction of a single crys- 
tal of TICUCI3 . Using a capacitance dilatometer we stud- 
ied the length changes perpendicular to the (010) and 
(102) cleavage planes of the monoclinic crystal struc- 
ture |5j . Via thermodynamic relations we derive the uni- 
axial pressure dependencies of the transition tempera- 
tures T c , of the transition fields H c , of the spin gap 
A, and of the magnetic coupling constants. For A we 
find huge pressure dependencies of about ±185%/GPa 
for uniaxial pressure perpendicular to the (010) or (102) 
planes, respectively. The uniaxial pressure dependencies 
of A unambiguously correlate with changes of the in- 
tradimer coupling J under pressure. In contrast to recent 
assumptions [lfj , pressure-dependent changes of the in- 
terdimer coupling J' play a minor role. Thus, the weaker 
J' of KCUCI3 is not a consequence of chemical pressure. 
This gives clear evidence that the Tl + and K + ions are di- 
rectly involved in the superexchange which is responsible 
for the relevant interdimer coupling and experimentally 
confirms the theoretical result of a significantly stronger 
superexchange via Tl + than via K + for i?CuCi3 . 

In Fig. ^ we show the longitudinal thermal expan- 
sion a.i along different directions i for various values of 
a magnetic field applied parallel to the respective mea- 
surement directions. In low fields there are no anomalies 
of oti, whereas above 6 T strong anomalies of opposite 
signs appear for both directions. With increasing field 
the anomalies increase and systematically shift to higher 
temperature. These anomalies signal spontaneous distor- 
tions below the phase transition: a spontaneous elonga- 
tion perpendicular to the (010) and a spontaneous con- 
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FIG. 1: Thermal expansion an perpendicular to the (010) 
and (102) planes for different magnetic-field strengths applied 
parallel to the respective measurement direction. 



traction perpendicular to the (102) plane. 

Fig. [21 displays representative magnetostriction mea- 
surements at different constant temperatures for both 
measurement directions. The relative length changes 
ei — ALi(H) I Li as a function of field for i — (010) and 
(102) are again of comparable size but of opposite sign. 
The phase transition causes a sharp kink in as a func- 
tion of field. With increasing temperature these kinks 
shift to higher fields and cannot be observed anymore in 
the studied field range above about 9K. For H > H c , 
changes essentially linear with field, whereas for smaller 
fields or for higher temperatures et is proportional to H 2 
(see Fig. A kink in is typical for a second-order 
phase transition, which should give rise to a jump-like 
anomaly in the field-derivative dei/dH. As mentioned 
above, there is some indication for a first-order contribu- 
tion to the phase transition in TICUCI3 . A small region of 
coexisting phases around H c as proposed from the NMR 
data [l5T ] can be neither confirmed nor ruled out by our 
measurements of the macroscopic length changes. We 
can, however, exclude that there is a significant hystere- 
sis of about 0.5 T between the H c values obtained with 
increasing and decreasing ma gnet ic field as has been ob- 
served in an ultrasound study [lj] . In Fig.|2t we compare 
dei/dH obtained with increasing and decreasing mag- 
netic field. As expected for a second-order phase transi- 
tion there is a (broadened) jump at H c and both curves 
agree well with each other over the entire field range. 
Thus, any hysteresis of H c is restricted to less than our 
field resolution of about 50 mT. 

In Fig. we show the phase diagram obtained from 
our data (circles) together with a power-law fit (solid 
line) of the form (g/2)[H c {T) - H c (0)] oc T*. This fit 
yields $ = 2.6 and H c (0) = 5.6 T. The exponent is 
larger than the value of 2.1 obtained in Ref. [2fT| (for 
T < 4K), but agrees well with the result of Quan- 
tum Monte Carlo (QMC) simulations [2l| . According to 
a more recent QMC study 22], <I> sensitively depends 
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FIG. 2: Magnetostriction Ei = ALi/Li perpendicular to the 
(010) and the (102) planes at selected temperatures. Panel 
(a) shows tiVS.H and (b) tiVS.H 2 , respectively. Panel (c) 
displays the field-derivatives for i = (010) obtained with 
increasing (line) and decreasing (symbols) field at T = 3.8 K. 



on the temperature range of the fit and in the low- 
temperature limit $ = 1.5 is approached, which agrees to 
the expected value for a BEC. Thus, our larger value of 
<!> arises most probably from the used temperature range 
(3K< T < 9K), but one should also keep in mind that 
$ could change due to the finite spin-phonon coupling, 
which is not considered in the models 0, @, 0, H3| • 

The anomalies at the phase boundary allow to derive 
the uniaxial pressure dependencies of T c and H c by the 
Ehrenfest relations 
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Here, V m is the molar volume, Acti is the height of the 
thermal-expansion anomaly (see Fig. ^) and AC P that 
of the corresponding specific- heat anomaly [jfl] . is 
the slope change of ei at H c (see Fig.EJ) and A the 
corresponding slope change of the magnetization |^|. 

With C p and M from Refs. 0, HJ we find huge uni- 
axial pressure dependencies of T c and H c , e.g. g p ^ c 1Q — 
-9K/GPa for T c = 7.2 K and H = 12 T, or ~ 

-8T/GPa for H c = 6T and T = 4K. A hypotheti- 
cal uniaxial pressure of 0.1 GPa on (010) would strongly 
shift the phase boundary towards higher H c and lower T c 
values, respectively, whereas uniaxial pressure on (102) 
would cause a shift in the opposite direction as shown 
in Fig. by the upward and downward triangles, respec- 
tively. The dashed and dotted lines represent power-law 
fits (keeping $ fixed) of these hypothetical phase bound- 
aries under uniaxial pressure, and their extrapolations to 
T = K reveal the uniaxial pressure dependencies of A 
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FIG. 3: Phase diagram of TIC11CI3 from the magnetostriction 
and thermal expansion measurements perpendicular to the 
(010) (o) and (102) (•) planes normalized by the respective 
g factors |2(3| . The triangles display the hypothetical shift of 
selected T c and H c values for a uniaxial pressure of 0.1 GPa 
perpendicular to (010) (A) and (102) (?) calculated by Eq.Q] 
The lines are power-law fits of the phase boundaries. 



These are huge values, but due to the opposite signs, 
they almost cancel each other under hydrostatic pressure. 
Nevertheless, a strong decrease of A has been observed 
under hydrostatic pressure 0, 0] . Thus A should also 
strongly decrease for uniaxial pressure along the [201] 
direction, which is perpendicular to both directions of 
our measurements. The geometry of our crystal did not 
allow measurements along the [201] direction, but we ex- 
pect that there will be similar anomalies at the phase 
boundary as those perpendicular to the (102) plane. 

In a model of dimers coupled by an effective interdimer 
coupling J' the magnitude of A is determined by the bal- 
ance of J and J'|lO|. An increase of J will enlarge A, 
whereas an increase of J' will enhance the bandwidth 
of the triplet excitations and therefore lower A. Due to 
the small value of A compared to J (and J') in TICUCI3 
already moderate pressure-dependent changes of J (or 
J') may cause drastic changes of A. In order to gain 
information whether these changes arise from a pressure- 
dependence of J or of J', we fit the magnetic suscepti- 
bility x UHl f° r temperatures well above the gap by 
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Here, Xo(T) is the susceptibility of non- interacting spin 
dimers with intradimer coupling J, and xmf accounts 
for a mean-field correction with an effective interdimer 
coupling J'. As shown by the solid line in Fig. 0] the 
fit for T > 25 K yields a good description of the exper- 
imental data for J = 60 K, J' = 53 K, and g — 1.48. 
Our value of J is close to the neutron scattering result 

j ~ 64 k H03H EH, whereas our J' is significantly 



FIG. 4: Magnetic susceptibility (o) with a fit for coupled 
dimers (solid line, Eq.|5J. The dashed and dotted lines show 
the expected changes of \ for an increase of J and a decrease 
of J', respectively. Both cases would cause an increase of the 
spin gap A (see text.). 



larger than the largest interdimer coupling (37 K) and 
our g factor is significantly smaller than (7010 — 2.06 ob- 
tained by ESR 20]. One has to expect quantitative dis- 
crepancies due to our oversimplified model. However, 
this does hardly affect the following analysis of the rela- 
tive variations around the maximum of x(T) arising from 
pressure-dependent changes of J or J', because the main 
result is obtained from the signs of the uniaxial pressure 
dependencies of x an d A, respectively. 

Let us discuss the case that A increases as it would 
under uniaxial pressure on (010). This may result ei- 
ther from an increase of J or from a decrease of J'. For 
both cases x(T) can be modeled by Eq.0 As shown 
in Fig. ^ the maximum Xmax decreases if J increases 
(dashed line) whereas Xmax increases if J' decreases (dot- 
ted line). Thus, the sign of the uniaxial pressure depen- 
dence of Xmax allows an unambiguous decision whether 
the uniaxial pressure dependence of A results from a 
change of J or of J'. Measurements of x under uniax- 
ial pressure do not exist. However, the uniaxial pressure 
dependence of x is related to the magnetostriction by a 
Maxwell relation, and £j cx H 2 is expected for a param- 
agnetic material with M = xH, i.e. 
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As shown in Fig. the relation oc H 2 is indeed ful- 
filled and for the discussed case one finds e io > (upper 
panel), which implies dxmax/dpoio < 0. A decreasing 
Xmax means that the intradimer coupling J increases (see 
Fig-EJ). The same argumentation with inverted signs for 
all the uniaxial pressure dependencies is valid for pres- 
sure on (102). Thus, J is the relevant parameter which 
changes under pressure! 

The anisotropy of dx ^ ax for i = (010) and (102) agrees 
well with that of SA. This indicates that J' hardly 
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changes under pressure Q. Setting = 0, our model 
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The relative changes of J are much smaller than those of 
A, but because A is much smaller than J, the absolute 
changes of J and A are not too different (±22 K/GPa and 
±14K/GPa, respectively). This can be interpreted as 
follows. The pressure- induced change of J causes mainly 
a shift of the center of mass of the triplet excitations, 
but hardly changes its bandwidth. This is completely 
different from what is observed when Tl is substituted 
by K: The triplet excitations of KCUCI3 have a much 
smaller bandwidth than those of TICUCI3. Our analy- 
sis of the pressure dependencies clearly shows that this 
smaller bandwidth is not a consequence of chemical pres- 
sure, although K + is significantly smaller than Tl + . Thus 
the very different values of J' in KCUCI3 and TICUCI3 
mean that the K + and Tl + ions directly influence the 
effective interdimer coupling. This conclusion has been 
proposed also from a bandstructure calculation 0] and is 
now experimentally confirmed by our data. The smaller 
J' arises most probably from a weaker overlap via the 
small [Ar] shell of K + than via Tl + with the configura- 
tion [Xe]4/ 14 5d 10 6s 2 . 

Although it is, in general, difficult to predict the mi- 
croscopic changes under (uniaxial) pressure, one may un- 
derstand qualitatively the uniaxial pressure dependen- 
cies of J in a simple microscopic picture. The dimers 
are formed by the Cu 2+ spins of two neighboring CuClg 
octahedra, which are connected via a common edge of 
their basal planes. The Cu-Cl-Cu bond angle amounts 
to ~ 96°. Thus the weak antifcrromagnetic coupling of 
TICUCI3 agrees with the expectation of the Goodenough- 
Kanamori- Anderson rules that the coupling changes from 
weakly ferro- to strongly antiferromagnetic when the Cu- 
Cl-Cu bond angle increases from 90° to 180°. Since the 
line connecting the two Cl~ ions and the [010] direction 
have an angle of about 29°, one may expect that pressure 
along [010] will shorten the Cl-Cl distance. This would 
increase the Cu-Cl-Cu bond angle and enhance J. The 
opposite may be expected for pressure along [201], since 
this direction has an angle of about 25° with the connec- 
tion of the Cu 2+ ions, and a shortening of the Cu-Cu 
distance would lower the Cu-Cl-Cu bond angle and de- 
crease J. The normal of the (102) plane is nearly perpen- 
dicular to both, the Cu-Cu (~ 82°) and the Cl-Cl line 
(~ 77°). Thus one may expect that pressure on the (102) 
plane will hardly change the Cu-Cl-Cu bond angle, but 
will slightly increase both the Cu-Cu and the Cl-Cl dis- 
tance. Therefore the Cu-Cl distances will increase and J 
decreases, since the overlap between the Cu-3d and Cl-2p 
orbitals becomes weaker. 

In summary, we have presented high-resolution mea- 



surements of thermal expansion and magnetostriction 
perpendicular to the (010) and (102) planes of TICUCI3. 
For both directions the field-induced Neel order causes 
very pronounced anomalies, which allow a detailed de- 
termination of the phase boundary. There is essentially 
no hysteresis as expected for a second-order phase transi- 
tion. The anomalies signal huge uniaxial pressure depen- 
dencies of the phase boundary, e.g. ±185%/GPa for the 
spin gap obtained from H c (0K) with the signs depend- 
ing on the direction of pressure. Large uniaxial pressure 
dependencies of opposite signs are also present for the 
susceptibility around 40 K. Our analysis unambiguously 
reveals that the huge pressure dependencies of A arise 
from pressure-dependent changes of the intradimer cou- 
pling, whereas changes of the interdimer coupling play 
a minor role. Thus, the smaller interdimer coupling in 
KCuCls compared to TIC11CI3 is clearly not a conse- 
quence of chemical pressure. 
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Since x-max and A depend on two parameters J and J', 
the anisotropies of a? l"'° x and of 4^- will be identical in 

opi apt 

two cases: (i) if pressure changes only one parameter (in 
our case J) or (ii) if the anisotropies of J^- and of 
are identical. Since the second case would be accidental, 
it appears very unlikely. 
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